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Abstract. The Wigner-Weisskopf-type model describing the energy transfer between 
two centers mediated by a continuum of energy levels is studied. This work is motivated 
by the recent interest in transport phenomena at nanoscale in biology and quantum 
engineering. The analytical estimation for the energy transfer efficiency is derived 
in the weak coupling regime and the conditions for the almost perfect transfer are 
discussed. The embedding of the standard tight-binding model into Wigner-Weisskopf 
one which includes the environmental noise is presented. 
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1. Introduction 

The almost perfect energy transfer (APET) in biologically relevant molecular systems 
attracts attention of physicists in the recent years [1]. Similar phenomena should 
be observed also in the engineered systems relevant for information processing like 
arrays of quantum dots or Josephson junctions [2j. Most of the theoretical analysis 
is performed using tight-binding Hamiltonian for an interacting N-body system in 
the single excitation regime. The influence of environment is usually modeled by 
Markovian master equations with phenomenological decay parameters. The numerical 
computations show the interplay between quantum propagation, quantum localization, 
decoherence and dissipation [H [Sj [5]. The aim of this paper is to present a simple, 
essentially exactly solvable model which describes a generic class of such phenomena and 
allows to derive analytical bounds on the efficiency of energy transfer processes. The 
motivation for the applied formalism comes from the particular model of two identical 
atoms placed in the focuses of two parabolic mirrors separated by a distance i, see Fig.l. 
We assume that the dipole moments of both atoms are parallel to the symmetry axis of 
the system. If initially the first atom is excited and the second is in the ground state 




Figure 1. The APET for two atoms places in thie focuses of two parabolic mirrors. 
Radiating dipoles are parallel to the axis. 

we expect that roughly after time t ~ i/c + 0{t), where r is a mean life-time of the 
excited state, the energy quantum carried by a photon is almost perfectly transfered 
to the second atom. The simple description of such a system can be given in terms of 
the Wigner-Weisskopf (W-W) Hamiltonian (compare the case of a single mirror and a 
single atom [7]). The conditions which allows for APET are encoded entirely in the 
spectrum of the Hamiltonian which is determined by the details of the model. Treating 
W-W model as a generic one of a large class of energy transfer phenomena between two 
localized centers one can try to find a general bound for the efficiency of such process 
and determine the sufficient conditions for APET. 

2. Wigner-Weisskopf model 

In order to introduce a 2-state W-W model we consider a quantum system of two 2-level 
atoms interacting with an electromagnetic field which can be described by the standard 
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Hamiltonian in the rotating wave approximation 

H = —al + —al + I d'kuik) at(k)a(k) 
2 2, J 

+ (^t o-iSi) + cria^gi) + (yt a(fi'2) + a^^a^g-i). (1) 

The operators ai^g^ and a^((7) are smeared annihilation and creation operators for the 
electromagnetic field given by a^{g) = f d^kg(k)a^(k), (7j(k),j = 1,2 are suitable 
formfactors, and the quantum fields satisfy [a(k), a'^(k')] = 5(k — k') (we omit for 
simplicity the polarization degrees of freedom). The standard 2x2 Pauli matrices 
refer to 2-level atoms. In the case of identical atoms with relative position described by 
the vector r the formfactors differ by the relative phase L(]s.) = k • r 

^2(k) = e-^«^7i(k). (2) 

The crucial property of the Hamiltonian ([T]) is the commutation with the excitation 
number operator defined as 

Ne. = J d^ka\k)aik) + {a', + 1/2) + {a'^ + 1/2). (3) 
The total Hilbert space of the system possesses a tensor product structure 

Htot = J^ph ® (4) 

with the bosonic Fock space J^ph describing an electromagnetic field. As a consequence 
of the commutation [H, N^x] = the subspace corresponding to the eigenvalue 1 of N^x, 
which can be called single exciton Hilbert space is invariant under the evolution. It 
possesses a direct orthogonal sum structure, i.e. is spanned by the following vectors 

|1) = I t;l) ® |vac) ® I ;;2), 
|2) =\i;l)® |vac) ® | t;2), 

I/) = I I; 1) (g) a^(/)|vac) | J,; 2) for any wave — packet /. (5) 

Here, | t; j) and | j) denote excited and ground state of the j-th atom, respectively, 
and I vac) is the vacuum state of an electromagnetic field. This mathematical 
construction, which has been used frequently in quantum optics [6], will be called the 
2-state W-W model. 

The 2-state W-W model is suitable for the following generic physical situation. 
We restrict ourselves to physical systems which can be described in terms of the single- 
exciton Hilbert space Tiex- In this Hilbert space we choose two orthogonal vectors |1), |2) 
which are determined by the initial state preparation and the measurement procedures. 
Namely, the initial state ("donor") of the system at time to = is denoted by |1), while 
the measurement at time t > after preparation is a von Neumann projection on the 
state 1 2) ("acceptor"). 

The Hilbert space of the model system is decomposed into a direct sum (compare 
with (^) 

nex = C®L\Q)®C (6) 
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where the one- dimensional subspaces C are generated by the states |1), |2) and the 
Hilbert space L^(f2) is their orthogonal supplement. For the convenience we represent 
this Hilbert space as the space of wave packets f{k),g{k), k ^ Q denoted by |/), \g) with 
the scalar product 



{f\9)= / f{k)g{k)dk. (7) 



n 



One should stress that the notation of above is used for convenience only. One can 
always replace the continuous variable A; by a joint set of continuous variables u and 
discrete quantum numbers m such that \k) is replaced by \co,m) and 

{uj,m\uj' ,m') = 5{uj — uj')5mm' , and / dA; i— > / dco. (8) 

The 2-state W-W Hamiltonian, which can be seen as the restriction of ([T]) to the single- 
exciton space, reads 

H =Ho + V, 

Ho = LOi\l){l\+co2\2){2\ + J dku{k)\k){kl 

V =(|l)(^?i| + k)(l|) + (|2)(^2| + |^?2)(2|) (9) 

where gi,2 € L^(f2). All details of the model are hidden in the form of formfactors 5^1,2 and 
the spectral resolution of the part of the Hamiltonian denoted here by dk u{k)\k) {k\. 
The only assumptions used in the following calculations are the continuity of the 
spectrum {uj{k)} and the orthogonality {gi\g2) = 0. The later condition means that 
there is no direct cross-talking between the preparation and measurement procedures 
and simplifies the calculations. 



3. Efficiency of energy transfer 

The system begins its evolution in the excited state |1) and after time t can be found 
in the state |2) with the probability 

Vu{t) = \Ai2{t)W ^i2 = (2|e-*^*|l). (10) 

The APET holds if Ui ~ U2 and for a certain time to the transfer probability 7^12 (^o) — 1- 
To find the sufficient conditions for APET in terms of the Hamiltonian notice that 
the (complex) probability amplitude Auit) is a scalar product of two states, e~*^*/^|l) 
and e*^*/^|2) which evolve forward and backward in time, respectively. For t ^ t where 
r is the life-time of the excited states one expects that the following approximation 
holds 

where fi{k), f2{k) are certain wave packets from the Hilbert space component L^(f2) 
describing intermediate excitonic states. The estimation (fTTj) is equivalent to the 
existence of wave operators defined as 

W+ = lim e'"°^e-'"^ , W_ = lim e-'"''^e'"^ (12) 
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which holds under mild conditions on the formfactors gi^2 and exciton's dispersion 
relation u}{k), at least in the weak sense of convergence of matrix elements. Combining 
(ini) with (Mh one obtains 



f,ik) = {k\W+\l) , f,{k) = {k\W^\2) 



and finally 



Ai2it) 



dke-'^^''^\k\W+\l){k\W-\2). 



(13) 
(14) 



ken 



The advantage of the Wigner-Weisskopf model is the fact that the matrix elements 
(l|iy_|fc) and (2|W4.|/i;) can be exactly computed using Laplace transforms (see the 
next Section). 



4. The computation of matrix elements 

The basic tools used for the computation of the probability amplitude (fHll 
the identity 

^.iHt ^ ^-iHot _ ^ [e-^"i(*-^)|l)((7i| + \g^{t - s)){l\ 

Jo 

+ e-*'^^(*-^)|2)((72| + \92{t - s))(2|]e-*^^ 
where \gj(t)) = e~^^°^\gj), and the Laplace transform 



are: 



e-''f{t) dt. 



(15) 
(16) 



(17) 



Introducing the notation (j,j' = 1,2): 

and using the identity f|T5|) . the definition of the wave operators ([T2|) . the notation ([1 
and the Laplace transform we have 

f,ik) = {k\W+\l) = - ^[g^ik)S^^i-^uik)) + g2ik)S2ii-icuik)) 

f^ik) = {k\W.\2) = t [giik)S2i Hco{k)) + g2{k)S22 Hcoik))] . (1; 
Combining the identity f lTSj) with the definitions f lT7|) we obtain a series of equations 



Suit) 


Jo 


-*'^i(*-^)Fii(s) 


ds, 


S22{t) 


Jo 


-''^'^'-'^F22{S) 


ds. 


Suit) 


= - i [ e-^'^i^*- 
Jo 


'^F,2{s) ds, 




Fiiit) 


= - ^ [ Gu{t- 
Jo 


- 5)5*11(3) ds, 
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F22{t)= - I f G22{t-s)S22{s)ds, 

Jo 

Fi2{t) = - ^ [ [Gu{t - s)Su{s) + Guit - s)S22{s)] ds, (19) 
which can be converted into equations for Laplace transforms 
~Sii{z) = ^—[l-iF,,{z)l 

Z + VjOi 

~S22{Z)=^—[1-IF22{Z)1 
Z + UX)2 

~Sx2{z) = ^^F,2iz), 

Z + tLOi 

Fn{z)= -zGu{z)Su{z), 

F22{z)= -iG22{z)S22{z), 

Fi2{z) = - t[Gn{z)Su{z)Gu{z)S22{z)]. (20) 
The system of equations (!20|) can be solved yielding the following basic formulas 

Su{z) = [z + iui + Gn{z)]~^, 
5*22(2;) = [z + iu2 + G'ii(z)]~\ 

Si2{z) = - Gi2{z)[z + lUi + Gii{z)]-\z + tuj2 + Gii{z)]-\ (21) 

Combining fl^ ) with flTSl) and f[T^ one can obtain the final formula for the transition 
amplitude which can be used for numerical calculations. 

5. Markovian approximation 

In order to continue analytical analysis of the problem we consider the case of weak 
coupling or Markovian approximation. Treating V given by ([9]) as a small perturbation 
one can omit in fllSp the terms proportional to 821- The further approximation concerns 
Laplace transforms Gjj{z). The real part of Gjj{—ujj) is a standard Fermi Golden Rule 
approximation for the decay rate of the state \j) (j = 1, 2) given by 



7j = Iji^j) , Iji^) = 71 dk \gj{k)\^5{uj{k) - u), (22) 



'n 

while the imaginary part is a radiative correction to the bare frequencies Uj. In the 
following we denote by the same symbol Uj the physical, renormalized frequencies. 
Summarizing, the Markovian approximation means that 

S2ii-iujik)) ~ , Sjj{-iuj{k)) ~ [i{ujj - uj{k)) + 7^-]"^ , < tuj. (23) 

Therefore using Eqs. [T8|) and fl23|) one obtains 



and 



m = .^/-w^. . = ^ (24) 

u{kj - Ui + Z7i u{kj - UJ2~ «72 



A,2{t)=l rffce-«*^-^ — (25) 

'fcgn (u;(fc) - wi +27i)(w(fc) -W2 + ^72) 
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One can estimate the upper bound for Vi2(t) using (12^ to get 

Pl2(t) = |A2(t)P< ||/2ini/l|P. (26) 

The equahty 7[x^ + 'y'^]~^dx = tt imphes for example 

ii/if = / \fm\'dk= [ dk 



n 



TT 



' .-1-^ (27) 



what finally gives the bound (compare with [8]) 

Vi,{t)<(l-^)(l-^)<l. (28) 
In the weak coupling regime the upper bound is close to one. 

6. Conditions for APET 

In order to approach the bound fl28|) and achieve APET certain matching conditions 
implied directly by the formula ( l25l) must be satisfied. To present them in a transparent 
form we introduce certain additional assumptions in a convenient parametrization. 
There are satisfied by the two-mirror system discussed in the Introduction. 

Assume that the intermediate Hilbert space L^(fi) is spanned by the basis \u,m) 
and the corresponding part of the Hamiltonian reads 

ill = / duju\u,m){u,m\. (29) 

The formfactors possess the following structure 

\gi) = / dug{uj)\u,mo) , \g2) = / due~'^'^'^'>g{u)\u,mo). (30) 



where g{oo) is a "fiat" slowly varying function and L{u) accounts for the spatial 
separation of donor and acceptor similarly to ([2]). They satisfy approximative 
orthogonality condition 

duje'^^''^\g{cu)\^ -0. (31) 

Now the mechanism leading to APET can be illustrated. We have to put the resonance 
condition 

001=002 = ooo, (32) 
which by (!22ll implies also the equality of decay rates 

7i = 72 = 7- (33) 
Then we have the approximative expression 

- i / d.- l^^(:^^^^)e<(^M-..). (34) 

71 J - Wo)"' + 7"^ Vcu - Wo + «7^ 
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Further approximations valid for |a; — WqI < 7 

^U-UJo^\ ^ _gi2(a,-a,o)/7 ^ ^(^) ^ ^(^^) + L'{uJo){uJ - UJo) (35) 
\UJ — LOq + / 

lead to 

\Ai2{t)\^-\ Idu- ^^-^—^exjp{i[L'{ujo) + 2T-t]{uj-ujo)}\ (36) 

TV J [UJ - Wo) + 7 

where r = I/7. The transition probability is maximal and close to one if the oscillating 
factor in the integral of above vanishes what happens for the time 

to = L'iujo) + 2t (37) 

which can be called optimal transfer time. 



7. Transport in quantum networks 

The standard tight-binding model of energy transport in quantum networks which is 
applicable both to molecular systems and engineered ones consist of N 2-level " atoms" 
described by the Hilbert space C and the Hamiltonian written in terms of Pauli 
matrices 

N N 

HN^^ujk(Tl+ (^ki (Tt^y; + h.c.) (38) 

k=l k<l=2 

where {cuk} are energies of the sites and {hki; k < 1} are hopping amplitudes. Similarly 
to the example in the Section 2 the exciton number operator 

N 

Ne. = + 1/2) (39) 

k=l 

commutes with Hj^^. 

Again a single exciton A^-dimensional Hilbert space is invariant with respect 
to the dynamics and the corresponding restriction of the Hamiltonian if at is given by 
the N X N hermitian matrix H = [hki] with the diagonal elements hkk — ^k- We 
attribute to the donor site an index " 1" and to the acceptor site index " 2" . We assume 
that the direct hopping 1 2 is neghgible, i.e. hi2 — h2i — and the two vectors 
\gi) = [0, 0, /i3i, /i4i, /iivi] and \g2) = [0, 0, /i32, /i42, /ijV2] are orthogonal Under 
these conditions the Hamiltonian can be recast into the discrete version of the W-W 
Hamiltonian 

H = a;i|ei)(ei| + uj2\e2){e2\ + {\ei){gi\ + \e2){g2\ + h.c.) + Hi. (40) 

Here |ei) = [1,0, ...,0], {62) — [0, 1,...,0] and Hi is a submatrix of H with indices 
k,l — 3, 4, N. Hi can be written in its spectral decomposition form 

N 

Hi^J2^-\a){a\ (41) 

a=3 
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to get a full analogy with the continuous W-W Hamiltonian In the next step one 
should embed the discrete model into a continuous one. Physically, it means that the 
interaction with other (environmental) degrees of freedom is taken into account. This 
interaction transforms the eigenstates \a) into resonances with finite spectral widths. 
The Hilbert space spanned by {|a)} is replaced by L^(]R) and the Hamiltonian (|4T]) by 
the multiplication operator by u which can be formally written as J duu \u){u\. As a 
consequence the vectors {\gj)]j = 1,2} are replaced by the wave functions gj{uj) which 
provide continuous envelopes for the discrete values {gji^a) = ■^j{'^\9j),<^ = 3, 4, N}. 
The constants N'j are chosen to satisfy normalization condition 



The construction of a continuous envelope for absolute values is simple as one 

can use e.g. Lorentz profiles to smear the eigenvalues {cq,}. However, the envelope of 
the relative phases denoted by L{u) is more tricky and depends on the geometry of the 
system (compare with ([2])). 

Using the results of the Section 6 we can discuss the conditions which should be 
satisfied to achieve APET, at least in the weak coupling regime. There are two types of 
such conditions, fine-tuning and generic ones. 
The fine-tuning conditions are : 

a) the resonance condition ( l32i) . 

b) the equality (l33ll which means that the formfactor functions gj{uj),j = 1, 2 must cross 
at the resonance energy uq. 

The generic conditions are the following: 

c) the optimal transfer time fl37|) should be short enough what implies that the relative 
phase function L{u) must be smooth enough around the value uo, 

d) the widths of the resonances should be larger than the nearest neighbor energy spacing 
for {ea} in order to produce flat envelopes gj{uj) what means that the interaction with 
an environment should be strong enough. 

In the case of engineered networks the fine-tuning conditions imply a careful design while 
for biologically relevant systems one can imagine that the natural selection mechanism 
plays a crucial role. 

8. Concluding remarks 

The 2-state Wigner-Weisskopf model provides a versatile mathematical tool to study, 
within the single exciton approach, the transport properties in complex molecular 
systems or engineered networks relevant for quantum information processing, both in 
the weak coupling (Markovian) and strong coupling (non-Markovian) regimes. The 
presented version of this model can be modified to include exciton decay and presence 
of a sink by adding imaginary decay rates iT{k) , iT2 to energies u{k) and 0J2, 
respectively. The W-W Hamiltonian (Q possesses also entirely classical interpretation 
as a Hamiltonian of a certain quadratic system (two harmonic oscillators coupled to a 




(42) 
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continuum of harmonic ones). Therefore, the physical mechanisms leading to APET are 
coherent wave-like ones and not specifically quantum which could be attributed to any 
nontrivial manifestation of entanglement. 
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